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Abstract 

The  Adaptive  Hydraulics  (ADH)  model  is  an  adaptive  finite  element  method  to  simulate  three-dimensional  Navier- 
Stokes  flow,  unsaturated  and  saturated  groundwater  flow,  overland  flow,  and  two-  or  three-dimensional  shallow- water 
flow  and  transport.  In  the  shallow-water  flow  and  transport,  especially  involving  multispecies  transport,  the  water 
depth  (h),  the  product  of  water  depth  and  velocities  (uh  and  vh),  as  well  as  water  depth  and  chemical  concentration 
(he)  are  dependent  variables  of  fluid-motion  simulations  and  are  often  solved  at  various  times.  It  is  important  for 
the  numerical  model  to  predict  accurate  water  depth,  velocity  fields,  and  chemical  distribution,  as  well  as  conserve 
mass,  especially  for  water  quality  applications.  Solution  accuracy  depends  highly  on  mesh  resolution.  Adaptive  mesh 
refinement  (AMR),  particularly  the  h-refinement,  is  often  used  to  add  new  nodes  in  the  region  where  they  are  needed 
and  to  remove  others  where  they  are  no  longer  required  during  the  simulation.  The  AMR  is  proven  to  optimize  the 
performance  of  a  computed  solution.  However,  mass  with  gain  or  loss  can  occur  when  elements  are  merged  due 
to  removing  a  node  at  mesh  coarsening.  Therefore,  we  develop  and  implement  the  mass-conservative  unrefinement 
algorithm  to  ensure  the  mass  conserved  in  a  merged  element  in  which  a  node  has  been  removed. 

This  study  describes  the  use  of  the  Galerkin  finite  element  method  to  redistribute  mass  to  nodes  comprising  a 
merged  element.  The  algorithm  was  incorporated  into  the  ADH  code.  This  algorithm  minimizes  mass  error  during  the 
unrefinement  process  to  conserve  mass  during  the  simulation  for  two-dimensional  shallow-water  flow  and  transport. 
The  implementation  neither  significantly  increases  the  computational  time  nor  memory  usage.  The  simulation  was 
run  with  various  numbers  of  processors.  The  results  showed  good  scaling  of  solution  time  as  the  number  of  processors 
increases. 

Keywords :  Adaptive  mesh,  finite  element  method,  mass  conservative,  shallow- water  flow,  and  ADH  model 


1.  Introduction 

Accurate  solutions  are  often  sought  when  we  solve  partial  differential  equations  (PDE)  for  scientific  problems.  A 
variety  of  numerical  errors  (e.g.,  spurious  oscillation,  numerical  spreading,  grid  orientation,  and  peak  clipping/valley 
elevating)  [1,  2,  3]  often  exist  when  we  use  numerical  methods  to  discretize  the  PDE  system.  These  errors  can 
be  partially  mitigated  through  the  implementation  of  higher  order  finite  element  or  finite  difference  discretizations, 
for  example,  TVD  and  ENO  methods  [4,  5],  with  the  increase  of  spatial  and  temporal  resolution  throughout  the 
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entire  domain  to  avoid  oscillation.  But  uniformly  refining  the  computational  mesh  throughout  the  domain  is  way  too 
expensive.  An  adaptive  mesh  refinement  (AMR)  approach,  which  is  also  called  h-refinement,  can  be  more  effective  in 
its  use  of  computer  resources  for  large-scale  applications.  The  AMR  improves  the  mesh  resolution  in  regions  where 
the  numerical  error  is  large,  while  keeping  the  mesh  coarse  elsewhere.  However,  the  AMR  technique  is  more  complex 
to  implement  but  benefits  many  applications  exhibiting  discontinuties,  shock  waves,  or  phase  changes  [6],  etc.  A  wide 
variety  of  AMR  approaches  have  been  explored  in  the  literature  [7,  8,  9,  10,  11,  12],  including  serial  and  parallel 
approaches  [13,  14,  15]. 

Mass  conservation  is  the  primary  requirement  when  AMR  takes  place.  In  fact,  AMR  includes  two  parts:  mesh 
refinement  and  unrefinement.  During  the  refinement  step,  it  is  straightforward  to  add  these  nodes  and  elements  in  a 
manner  that  preserves  mass,  such  as  linear  interpolation  in  a  linear  element.  For  unrefinement,  the  nodes  are  deleted 
and  elements  are  merged  in  some  regions  where  a  fine  mesh  resolution  has  become  useless  for  transient  problems.  The 
unrefinement  step  is  necessary  to  avoid  memory  and  computation  overkill.  Without  special  care  of  transferring  data 
from  the  unrefined  node  to  its  parents,  the  result  often  brings  up  mass  conservation  issues.  Yamoah  [16]  investigated 
this  problem  and  proposed  three  algorithms  to  restore  mass  conservation  for  unrefinement  when  solving  1-D  variably 
saturated  groundwater  flow  using  MATLAB  code.  These  methods  are  based  on  either  a  weighted  average  approach 
or  the  L2  projection  approach.  He  concluded  that  the  mass  cannot  be  preserved,  and  the  mass  lost  may  be  significant 
for  a  long  simulation  if  a  mass-conservative  unrefinement  algorithm  is  not  used. 

We  extended  the  L2  projection  algorithm  onto  2-D  unstructured  meshes.  We  developed  the  mathematical  formu¬ 
lation  and  implemented  the  algorithm  in  the  ADaptive  Hydraulics  (ADH)  modeling  system — a  parallelized  numerical 
model  with  the  following  components:  saturated  and  unsaturated  groundwater  flow,  three-dimensional  Navier-Stokes 
flow,  two-  or  three-dimensional  shallow- water  flow,  and  the  associated  transport.  When  refinement  is  on  the  ADH,  the 
element  is  split  if  the  solution  error  on  the  element  exceeds  the  refinement  error  tolerance.  At  the  added  node  the  data 
such  as  water  head  or  chemical  concentration  are  linear  interpolation  in  the  element.  When  the  solution  error  on  the  el¬ 
ement  is  below  the  unrefinement  error  tolerance  then  the  elements  are  recombined.  The  node  is  simply  deleted  without 
special  care  of  transferring  data  from  the  unrefined  node  to  its  parents.  In  this  paper,  we  first  discuss  the  shallow- water 
flow  application  including  the  governing  equations  and  discretized  formulation,  numerical  methods,  and  numerical 
difficulties.  We  then  present  the  mathematical  formulation  of  the  mass -conservative  unrefinement  algorithm,  followed 
by  its  implementation.  Experimental  results  are  then  presented  to  demonstrate  the  mass  conservation  throughout  the 
entire  simulation.  Finally,  we  summarize  and  discuss  these  results  and  also  present  plans  for  future  works. 


2.  ADaptive  Hydraulics  (ADH) 


ADH  is  a  modular,  parallel,  adaptive  finite  element  model  for  one-,  two-,  and  three-dimensional  flow  and  transport. 
ADH  simulates  saturated-unsaturated  groundwater  flow,  internal  flow,  and  open-channel  flow.  The  groundwater  mod¬ 
ule  includes  solving  3-D  variably  saturated  groundwater  flow  and  constituent  transport.  The  groundwater  flow  may  be 
coupled  to  2-D  diffusive  wave  or  dynamic  wave  surface  water  flow  and  can  be  loosely  coupled  to  3-D  subsurface  heat 
transport.  The  internal  flow  module  is  used  to  solve  3-D  Navier-Stokes  equations.  This  module  is  nonhydrostatic  and 
includes  a  k  -  e  turbulence  model.  There  are  two  approaches  for  open-channel  flows:  nonhydro  static  and  hydrostatic 
pressure  methods.  The  nonhydrostatic  approach  is  more  appropriate  for  domains  near  structures  where  the  vertical 
inertial  terms  are  significant.  This  approach  is  the  same  as  that  used  for  the  3-D  internal  flow,  though  the  free  surface 
is  updated  dynamically.  The  hydrostatic  approach  is  to  solve  the  shallow-water  flow  equations,  including  2-D  and 
3-D  variable  density  flow  (baroclinic).  The  2-D  module  includes  the  capability  for  wetting-drying,  dam-break,  super¬ 
critical  and  subcritical  flow,  sediment  transport,  and  bed  change,  and  also  has  the  capability  to  include  the  long- wave 
vessel  effects. 

The  2-D  shallow- water  flow  equations  are  obtained  by  vertically  integrating  the  mass  and  momentum  under  the 
assumptions  of  incompressible  flow  and  hydrostatic  pressure  [17].  Assuming  negligible  free  surface  shear  and  negli¬ 
gible  fluid  pressure  at  the  free  surface,  the  2-D  shallow-water  equations  are  written  as  the  following  partial  differential 
equation  (PDE). 


dQ  <^Fx 
dt  dx 


+  -t—  +H  -  0 

dy 


(1) 


where  Q  =  \h,  uh,  vh}T,  t  is  time,  Fx  =  {uh,  u2h  + 


gh 2 


-  h^,uvh  -  h^] 
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H  =  {0,  gh^  +  ghn2u^f3+—,  gh^  +  gh n2y^3+— }r,  p  is  the  fluid  density,  u  and  v  are  the  flow  velocity  in  the  x-  and 
y-direction,  respectively,  h  is  the  water  depth,  g  is  the  gravitational  acceleration,  Zb  is  the  riverbed  elevation,  n  is  the 
Manning’s  roughness  coefficient,  and  <t*7  is  the  Reynolds  stress  due  to  turbulence,  where  the  first  subscript  indicates 
the  direction  and  the  second,  the  face  on  which  the  stress  acts.  The  Reynolds  stresses  are  determined  by  the  gradient 
in  the  mean  currents  using  the  Boussinesq  approach 


du 


dx 


dv 


dy 


=  2Pn  hr  *  °~yy  =  2Pyt  \- r)^xy  =  cr  =  2 pvt  —  +  T"  ’ 


du  dv 


dy  dx 


(2) 


where  vt  is  the  kinematic  eddy  viscosity,  which  varies  spatially. 

The  PDE  system  in  (1)  can  be  solved  using  the  finite  element  method  with  the  approach  of  Petrov-Galerkin  that 
incorporates  a  combination  of  Galerkin  test  function  and  a  non-Galerkin  component  to  control  oscillations  caused  by 
convection  [18].  The  variational  forms  of  (1)  in  integral  form  and  in  finite  element  form  are  given  as  follows: 


X 


SQ  +  3F* 
dt  dx 


dFy 

dy 


+  H 


d£l  =  0, 


dt  dx 


=  0, 


(3) 

(4) 


where  Q  is  the  domain,  the  subscript  e  indicates  a  particular  element,  the  subscript  i  indicates  a  particular  test  function, 
and  vt/i  is  the  test  function.  The  weak  form  of  the  equation  is 


Z 


Fy 


<PiA§B§  +  'F,H)<«2e+ 

J|.  [<Pi  (F xnx  +  Fyny)]  dTe 


=  0, 


(5) 


where  is  the  Galerkin  test  function,  tpt  is  the  non-Galerkin  part  of  test  function,  i.e.,  fi  ■  [Ax^fA  +  Ay^Z?],  Ax  and 
Ay  are  the  grid  interval  in  the  x-  and  y-direction,  respectively,  A  =  B  =  ,  re  is  the  boundary,  and  nx  and  ny  are, 

respectively,  the  x-  and  y-component  of  the  unit  outward  vector  normal  to  the  boundary.  A  description  of  A  and  B  can 
be  found  in  [19].  The  temporal  derivatives  are  expressed  using  a  finite  difference  expression  as 


\  dt  )  j  5tk 


a -«>•[***],• 


(6) 


where  the  subscript  j  is  the  node  index,  k  is  the  time-step  index,  and  a  is  the  temporal  relaxation  factor.  A  complete 
description  of  the  time-stepping  methods  employed  in  ADH  can  be  found  in  [20]. 


3.  Mathematical  Formulation  for  Mass- Conservative  Unrefined  Element 


In  this  section,  we  present  the  mathematical  formulation  of  the  mass-conservative  unrefinement  algorithm.  The 
derivation  is  based  on  the  backward  Euler  time  scheme.  The  time  derivative  term  in  (4)  can  be  formulated  using  the 
Galerkin  finite  element  method  based  on  a  set  of  linear  basis  functions  spanning  a  finite  dimensional  space.  Selecting 

N  N 

a  test  function  'F  =  2  Pi^i  and  finding  a  set  of  qt  for  the  trial  function,  we  can  approximate  Q  =  2  97 Qi  and  then 

i=  1  i= 1 


make 


For  AMR  problems,  the  mesh  may  vary  from  time  to  time  resulting  in  the  set  of  basis  functions  qn  e  £ln  and 
qn+l  e  Qw+1  at  time  tn  and  tn+l ,  respectively.  The  associated  finite  dimensional  space  functions  are  denoted  by  V£  and 
V"+1  at  two  consecutive  time-steps.  Given  Q"  e  V",  we  need  to  find  Q”+1  6  V"+1  such  that  the  backward  Euler  time 
difference  of  (7)  becomes 
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Qhn+lx¥hn+lda-  f 

jQn+1 


Q  hV+ldQ  =  0 


for  all  fVi+l  6  V/,"+i .  The  second  term 


(8) 

(9) 


introduces  problems.  The  function  Qhn  is  defined  on  the  mesh  at  previous  time  tn ,  but  the  test  function  is  defined  on 
the  domain  £ln+l  at  time  tn+l .  In  the  case  of  unrefinement,  Vhn+l  c  Vhn  and  Q ^  cannot  be  expressed  on  Qw+1  exactly. 


3.1.  Mass -conservative  merged  element  using  local  approach 

The  goal  of  the  mass  conservation  scheme  is  to  preserve  the  mass  in  a  merged  element,  in  which  a  node  has 
been  removed.  In  Figure  1,  elements  e\  and  e 2  are  merged  to  element  e*  due  to  unrefining  the  node  (*4, 34).  The 
local  approach  to  accurately  calculate  (9)  allows  a  temporary  and  possibly  discontinuous  solution  Q*  for  the  merged 
element  only.  Therefore,  compute  the  {QJ,  QJ,  Q3}  values  using  the  basis  functions  {TV}  such  that 

Le,  (Qi-W*1  +  Q2*N2n+1  +  Q3*N3n+x)%n+xdA  = 

Lei  (Qi "M"  +  QVW  +  Q3nN3n)'¥"+1dA  + 

L  (Q4nNin +  Q2nN2n +  Q3nN3n)'Vin+1dA  (10) 

JlLe2 


Figure  1 :  Element  with  a  unrefined  node 


3.2.  Local  linear  system  to  solve 
3.2.1.  LHS 

On  the  left-hand  side  of  (10),  the  formulation  can  be  rewritten  as  shown  in  (1 1),  and  the  integration  is  as  shown  as 
(12),  where  Ae*  is  the  area  of  element  e*. 
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Nin+1  ' 
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LHS  = 
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q*2 

Jne* 
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1  J  1 

{  q;  J 

(id 


LHS  =  T§ 


2  1  1 
1  2  1 
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Q; ) 

Q2  •  (12) 

q;  J 


3.2.2.  RHS 

On  the  right-hand  side  of  (10),  the  formulation  can  be  rewritten  to 
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/ 

where  (Q)  =  [Q”  Q"  Q'(  Q"|r. 

Once  the  mathematical  form  of  (13)  is  obtained,  the  solution  of  {Q*}  can  then  be  solved  using  a  direct  solver  from 
Equations  (12)  and  (13). 
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4.  Software  implementation 

The  data  structure  is  designed  as  in  the  following  box  on  the  left  to  store  the  merged  element  number,  the  edge 
number  that  the  unrefined  node  sits  on,  the  unrefined  nodal  coordinate,  the  quantities  stored  on  the  unrefined  node, 
and  the  {Q*}  values  solved  by  the  UNR  algorithm  described  in  the  previous  section.  Each  merged  element  needs  152 
bytes  of  memory  to  store  information  for  computation  in  order  to  conserve  water  mass  for  flow  simulations. 

The  following  box  on  the  right  depicts  the  implementation  of  solving  the  shallow-water  equations  with  AMR 
and  mass -conservative  unrefinement  algorithm  implementation.  The  function  adpt  jnain  is  the  main  driver  for  AMR 
comprising  (1)  mesh  unrefinement  (adpt_unref),  (2)  sweeps  of  mesh  refinement  (adpt_ref),  and  (3)  repartition 
if  needed.  At  the  end  of  each  of  these  two  functions,  node_renumber  and  elem_renumber  are  called  to  renum¬ 
ber  nodes  and  elements.  The  UNREF_elem  data  structure  starts  to  fill  in  for  each  merged  element  in  the  function 
adpt_unref ,  though  ustar  is  computed  later  in  the  function  UNRcompute.  Note  that  element  merge  is  not  allowed 
if  this  is  a  partial  wet  and  dry  element.  If  the  merged  element  is  split  again  in  the  function  adpt_ref ,  the  data  struc¬ 
ture  for  this  element  needs  to  deactivate  by  setting  the  field  ielem  to  -1.  The  data  structure  also  needs  to  keep  up 
with  the  new  element  number  assigned  in  elem_renumber.  If  the  merged  element  is  sent  to  another  processor,  the 
UNREF_elem  data  must  be  packed  for  the  send  and  unpacked  after  the  receive  on  the  receiving  processor.  The  algo¬ 
rithm  for  the  mass-conservative  unrefined  element  described  in  the  previous  section  is  implemented  in  UNRcompute. 


5.  Experimental  Results 

Two  ADH  models,  with  and  without  a  mass-conservative  unrefinement  algorithm  (UNR),  were  built  for  investi¬ 
gating  mass  errors.  The  results  of  the  volume  mass  obtained  from  running  applications  using  these  two  models  are 
compared  at  each  simulation  time-step.  Two  2-D  test  cases  are  selected:  (1)  tidal  flow  with  only  water  flow  and  (2)  a 
closed  pond  system  with  water  flow  and  conservative  constituent  transport. 

5.7.  Tidal  flow  test 

Figure  2  shows  the  initial  unstructured  mesh  containing  1843  nodes,  each  with  three  degrees  of  freedom  described 
in  (1),  and  3424  triangular  elements.  Water  flowed  in  at  the  top  and  out  through  the  left.  The  right  side  of  the  bottom 
boundary  is  a  dead  end.  The  time- step  size  is  600  seconds  throughout  the  entire  simulation  for  240  hours.  The 
convergence  tolerance  for  nonlinear  iterations  is  set  to  10-6,  and  the  maximum  number  of  nonlinear  iterations  is  set 
to  10. 

We  first  run  the  application  with  no  refinement.  The  initial  volume  mass  is  8.560000000000000  x  108  m3. 
The  water  mass  error  is  almost  close  to  zero  at  each  simulation  time-step  as  expected.  The  maximum  error  is 
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Figure  2:  Tidal  flow  test  domain  and  its  discretized  unstructured  mesh 


5.125999450683594  x  10-6  m3,  which  is  close  to  the  machine  error,  i.e.,  the  order  of  10-6  out  of  hundreds  mil¬ 
lion.  We  then  run  the  same  case  but  with  adaption.  The  refinement  tolerance  (f)  is  1  m3/s,  and  the  unrefinement 
tolerance  is  0.001  m3/s.  At  each  time-step,  the  model  performs  adaptive  refinement  if  the  error  indicator  is  above  1.0 
m3/s  and  unrefinement  if  the  error  indicator  is  below  0.001m3/s.  Figures  3(a)  and  3(b)  show  the  number  of  merged 
elements  and  number  of  total  elements  at  each  time-step.  They  are  changed  only  at  the  early  and  late  stages  of  sim¬ 
ulation.  Figure  4  plots  water  mass  error  from  the  runs  with  or  without  the  UNR  implementation.  The  mass  error 
without  the  UNR  algorithm  ranges  from  -48.2196  m3  to  101.2599  m3,  while  with  the  UNR  algorithm,  it  ranges  from 
-7.1525  x  10-6  m3  to  8.9406  x  10-6  m3.  From  the  result,  the  UNR  algorithm  does  sustain  the  same  order  of  mass 
error  from  the  no-adaption  case.  We  can  also  observe  that  the  more  merged  the  elements,  the  worse  the  mass  error  is 
if  UNR  is  not  implemented.  That  also  explains  why  the  same  mass  error  is  maintained  as  the  simulation  approaches 
to  the  end.  Note  that  the  setup  of  this  test  is  an  easy  physics  problem  to  solve  so  that  the  mass  error  is  small  even 
without  mesh  adaption. 

5.2.  Pond  test 

The  pond  shown  in  Fig.  5  contains  1408  nodes,  each  with  three  degrees  of  freedom  described  in  (1)  and  2664 
triangles  initially.  No  inflow  and  outflow  occur  across  the  closed  system.  In  addition  to  the  flow  simulation,  one 
chemical  constituent  is  involved  in  the  transport  process.  The  time-step  size  is  10  seconds  for  the  entire  simulation 
of  6000  seconds.  The  convergence  tolerance  for  nonlinear  iterations  is  10-8,  and  the  maximum  number  of  nonlinear 
iterations  is  set  to  10. 

The  initial  water  volume  is  5.4298776489890  x  106  m3.  With  AMR  turned  off,  Figures  6(a)  and  6(b)  display  the 
water  mass  error  and  the  chemical  mass  error  vs  simulation  time.  The  water  mass  error  ranges  from  -2.14204192  x 
10-8  m3  to  2.51457095  x  10-8  m3,  while  the  chemical  mass  error  ranges  from  -3.35276136X  10-8  m3  to  3.35276126X 
10-8  m3.  With  mesh  adaption  turned  on,  the  simulations  are  run  with  and  without  the  UNR  algorithm.  For  both  runs, 
the  refinement  tolerance  is  1  m3/s,  and  unrefinement  tolerance  is  0.001  m3/s.  Figure  7  shows  the  number  of  total 
elements  and  the  number  of  merged  elements.  During  the  early  simulation,  the  total  number  of  elements  increases 
quickly  to  the  maximum  and  then  gradually  reduces  to  the  initial  value,  while  a  large  amount  of  merged  elements  can 
be  found  between  1000  and  3000  seconds  of  simulation  time. 
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Figure  3:  (a)  and  (b) 


Water  Mass  Err  init.  mass  =856000000Max  err=1 01 .2599 
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Figure  4:  Water  mass  error  with  and  without  the  UNR  algorithm 


Figures  8(a)  and  8(b)  show  water  mass  error  and  chemical  mass  error  from  the  run  using  AMR  with  and  without 
the  UNR  algorithm.  The  water  mass  error  without  the  UNR  algorithm  ranges  from  -14.0923  m3  to  0.2799  m3,  while 
with  the  UNR  algorithm  it  ranges  from  -3.4458  x  10-8  m3  to  3.1664  x  10~8  m3.  For  the  chemical  constituent,  the 
chemical  mass  error  without  the  UNR  algoritm  ranges  -177.8483  g  to  708.6960  g,  while  it  ranges  from  -4.0978  x  10-8 
g  to  2.2383  x  10-8  g  with  the  UNR  algorithm.  Therefore,  the  UNR  algorithm  can  reduce  the  water  and  chemical  mass 
error  to  the  same  order  of  magnitude  as  that  obtained  by  running  the  case  without  adaptation.  Similar  to  the  previous 
test  case,  the  pond  case  is  an  easy  system  to  solve  because  the  water  mass  error  is  within  the  range  of  the  machine 
truncation  error  even  without  mesh  adaption.  The  UNR  algorithm  is  used  to  maintain  the  mass  conservation  for  AMR 
cases,  which  cause  mass  loss  or  mass  gain  if  just  simply  deleting  a  node  during  unrefinement. 

Table  1  shows  the  total  wall  clock  time  in  seconds  vs.  the  number  of  processors.  The  application  was  run  with 
1,  2,  4,  and  8  processors  on  the  Cray  XT4  system.  The  application  only  ran  with  less  than  8  processors  because  of 
the  limit  of  the  small  total  number  of  elements.  The  result  depicts  the  overhead  of  using  the  UNR  algorithm  is  small. 
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Figure  5:  Domain  of  the  pond  and  its  initial  discretized  unstructured  mesh 
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(b)  Chemical  mass  error  without  adaptation 


Figure  6:  (a)  and  (b) 


Some  computational  operations  are  executed  in  the  UNR  routine  and  increase  the  total  wall  clock  time  compared  with 
the  case  without  the  UNR  algorithm.  However,  the  overhead  resulting  from  the  UNR  algorithm  is  small  and  can  be 
offset  by  the  reduced  number  of  nonlinear  iterations.  This  is  achieved  by  computing  more  accurate  solutions  each 
time. 


Table  1 :  Total  wall  clock  time 


number  of  processors 

Total  wall  clock  time  (secs) 

with  UNR 

without  UNR 

1 

1185.4990 

1175.0420 

2 

677.4126 

676.9971 

4 

351.0934 

352.3615 

8 

177.9629 

173.1149 

Hung  V.  Nguyen  et  al.  /  Procedia  Computer  Science  9  (2012)  727  -  736 


735 


Figure  7:  Total  number  of  elements  at  each  time  from  the  pond  case 
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(a)  Water  mass  error  with  and  without  the  UNR  algorithm  (b)  Chemical  mass  error  with  and  without  the  UNR  algorithm 


Figure  8:  (a)  and  (b) 


6.  Conclusions  and  future  works 

The  UNR  algorithm  is  developed  and  implemented  successfully  in  the  ADH  model  for  two-dimensional  shallow- 
water  applications.  It  is  known  that  there  is  mass  error  when  merging  elements  in  the  common  implementation  of  the 
AMR  procedure.  We  demonstrate  that  the  UNR  algorithm  can  reduce  water  and  chemical  mass  error  to  the  same  order 
of  magnitude  as  the  mass  error  obtained  by  running  without  mesh  adaptation.  In  addition,  the  overhead  incurred  in 
the  UNR  algorithm  is  considered  small.  The  implementation  requires  a  small  amount  of  memory  to  store  and  transfer 
information  on  merged  elements  when  running  on  multiple  processors.  This  research  currently  is  presented  for  the 
backward  Euler  time  scheme.  The  mass  error  is  actually  a  good  indicator  for  the  setup  of  mesh  resolution,  nonlinear 
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tolerance,  and  other  parameters  in  order  to  obtain  accurate  solutions. 

We  are  planning  to  implement  the  UNR  algorithm  with  the  higher  order  time  scheme  in  the  ADH  model  for  two- 
dimensional  shallow- water  flow.  Although  all  the  demonstration  cases  are  small-scale  problems,  the  extension  to  large 
scale  is  without  limits  and  only  requires  data  set  generation.  The  application  could  also  be  extended  to  other  modules 
available  in  ADH  such  as  3-D  variably  saturated  groundwater  flow.  Mathematical  formulation  for  3-D  applications 
will  soon  be  discussed. 
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